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ABSTRACT 


A light scattering apparatus to measure particle size 


| Solid rocket motor was improved. Multiple 
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consecutive scans of two photodiode arrays were accompiished 


with a pacing circuit and added memory. The device was 


calibrated using various suspended particle samples and 


found to make accurate measurements. 
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[аташе се predic tIon codes for solid rocket motors 
model two phase flow losses as functions of particle size. 
ijeadceition, particle size within the grain Sort Strongly 
atfects the damping of combustion pressure oscillations. At 
present these models are based on particle size data from 
collected exhaust samples [Ref. 1]. However, particle size 
varies with position in tne motor and otner parameters 
(pressure, propellant formulation, nozzle design, etc.). 
Therefore, experiments to determine how particle size varies 
Inh actual flow environment of the motor (i.e., across the 
nozzle) are needed to validate the models for two phase flow 
losses. Cramer [Ref. 2] and Karagounis [Ref. 3] provide a 
good summary of the subject and the Naval Postgraduate School 
Memoustion Laboratory effort to obtain particle size data 
across the exhaust nozzle of a solid rocket motor. 

Еее Ева та this continuing e fort was the dif- 
fractively scattered light technique. The diffraction patterns 
of light scattered by particles are analyzed to determine the 
volume to surface mean diameter [Refs. 4 through 11]. This 
method has the disadvantages that size distributions cannot 
be easily determined and particles larger than some threshold 
Ese will not be Бе. to the ЕЕЕ шау small angles 


as which they scatter light. However, it has the advantage 


that it is non-iatrusxus апсы сес = К - c 
internal motor environment: 

Use of this method for particle sizing in solid roc и 
motors at the Naval Postgraduate School Combustion Labord ea 
was begun by Karagounis [Ref. 3l. The apparatus was sub- 
sequently redesigned and the data acquisition equioment up- 
graded with the introduction oz the Hewlett Packard 3054A 
data acquisition system with an HP 85 as the controlling 
computer [Ref. 12]. The investigation by Cramer and Hansen 
followed and showed that propellent composition can limit 
the application of the technique. Large particulate combusibeem 
products in the flow made particle size data difficult to 
obtain. This was especially true if only one measurement of 
the scattering profile was made during a test firing. 

To address this problem in the present study several 
improvements were made. A cleaner burning propellant was 
obtained to reduce char agglomerates in the exhaust producte 
A more statistically valid data sample (multiple measurements 
during a single test) was made possible with added momory 
in the data acquisition equipment and a pacing circuit which 
allowed full use of this memory. Data reduction was also 
improved with a Hewlett Packard 9836S computer combined with 
a more recently developed approach to particle sizing 
presented by Buchele [Ref. 13]. е method is discussed 


later in depth. 
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Implement the HP $8365 as the system controller. 

Expand the multiprogrammer memory in order to obtain up 
to eight consecutive scans of the diode arrays during 

A4 EeSt ۲3 210 ٠ 

Improve data reduction techniques by the method of 
Beehele (Ref. I3]. 

Modify the apparatus and experimental procedures to 
improve the angular resolution and to reduce extraneous 
Диска 

Certify the proper functioning and accuracy of the 


s ЛО DEIOI tO actual motor testing. 


TI 


ТІ. ۷٢٣ Е 


А, GENERAL Discuss | 

The completely general theory of scattering was developed 
by Mie and is presented by Van de Hulst [Ref. 14]. Тһе 
light scattering characteristics for spherical particles of 
any size are fully described in a mathematical format. The 
Mie scattering functions contain Legendre polynomials and 
spherical Bessel functions and fully treat the phenomena of 
reflection, refraction, diffraction, and Sexfinceicn. sua. 
full theory is most often applied when particle size is 
approximately the same as the wavelength of the incident light. 
Van de Hulst [Ref. 14] calls this the regime of Anomalous 
басе 7۳۴ 

For particle sizes much smaller than the wavelength of 
light the Mie equations simplify to a form which is more 
dependent on the index of refraction of the particles and 
less dependent on particle size. This is called Rayleigh 
Scattering. 

The study of particle size behavior in solid propellant 
rockets mainly covers sizes much greater than the wavelength 
of light. Scattering by large particles such as these is 


described adequately by Fraunhofer diffraction. 
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The ringed diffraction pattern aenerated by a hole in 
a mask, or a number of particles of the same size is de- 
scribed by the equation: 


2 
2J4 (a6) 


| 
a8 


where: 
I(0) describes the relative intensity of the scattered 
light at an angle theta (9) 
J4 (a9) is the Bessel function of the first kind 


TD А : ۱ 
ф = = is the particle size parameter for diameter 


À 

D and wavelength of light lambda (A). 

Measuring the particle size for a monodispersion can be 
accomplished by measuring the angular position of a dark or 
Meme rine in the diffraction pattern. For a dark ring the 
zero of the Bessel function corresponding to the ring is set 
спа расе Size teudetemmaned directly [Ref. 5]. 
For bright rings one sets (a0) equal to the corresponding 
maximum of the Bessel function and solves for the diameter. 

The above method is not used for polydispersions since 
the discrete rings are not observed. However, Dobbins, et al. 
[Ref. 5] introduced a significant improvement in the dif- 


асе ме scattered light method of particle sizing. They 


found that although the method was not directly able to 
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mean diameter defined by 


pan (г) пар 

P Б 

Dae = 5 (1) 
2 *N. (D)D dD 


су С) 


where: 

N (D) is a distribution function describing the pro- 

portion of particles with diameter (D) in the sample, 

could be accurately measurec. 
A curve for sizing polydispersions was presented which was 
used by Cramer and Hansen [Refs. 2, 12]. 

Two phase flow losses are often calculated in terms of 
043: If the distribution of sizes in the polydispersion is 
well behaved then 042 and other diameters such as D 13 can 
be easily related [Ref. 15]. Reference 5 reported that very 
small particles in the distribution have a minor influence 
on the scattering profile. This makes the measurement 
technique promising for the two phase flow loss study 
since very small particles do not contribute significantly 
to these losses and so are of less interest. 

Roberts and Webb [Ref. 6] essentially confirmed the 
conclusions of Reference 5 and presented a similar curve for 
use in sizing. 

More recently, Buchele [Ref. 13] gives a good summary of 


experimental techniques for particle sizing by measuring 
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КЕ он о и са: 52205505 С: 1255255 12 2165 recort 1 
шшаси ше тергезепъз the seasctering prorile of a polvdispersion 
with a function which closely approximates the curves of 


References 5 and 6. 


I (8) = EXP -(.5708)* 


This function from Reference 13 and the curve from 
Reference 5 were both used in the present study to evaluate 
the apparatus to be used with solid propellant rocket motors. 

An additional detail of measuring scattering profiles is 
covered by Van de Hulst [Ref. 14]. The wavelength of light 
used in the scattering calculations depends on the index of 
refraction of the medium containing the particles. The 


wavelength used in all calculations must be: 


met 
M 
where: 
Ag is the wavelength of light in a vacuum and, 


М is the refractive index of the medium with respect 
O a vacuum. 


Thus, the size parameter (a) becomes: 


TDM 


Ло 





Yo = 


and the beam spread parameter becomes: 
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POM 


2 


(со) = 


Another consideration is presented by Gumprecht and 
Sliepcevich [Ref. 4]. Light scattered by particles in a 
medium is refracted as it crosses each interface of the 
container holding the medium. This is discussed further in 
the section on calibration and evaluation of the apparatus. 

Additional complications arise with the full treatment 
of the index of refraction of the particles with respect to 
the medium. But, for Fraunnoter diffraction alone this aspect 


can be neglected. 


C. RESTRICTIONS AND SOuRGChsw of Sennen. 

Some restrictions on use of the method are described by 
Dobbins, et al. [Ref. 5] and were satisfied as described by 
Cramer [Ref. 2]. These are related to the size of particles, 
the distance to the detector, and some phenomenon covered in 
the rigorous Mie theory. 

One must keep in mind also that the curves developed for 
polydispersions are based on the Upper Limit Distribution 
Function of Mugele and Evans [Ref. 7]. This means that no 
particles with size greater than approximately ten times the 
mean should be in the sample [Ref. 13]. This appears to be 
a mild restriction. Van ае ние [ке um Ue 


criteria for single Scattering and a simple test to verify 
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prooortional to the number of particles the mathematics 
remain simple. 

Sources of error of the diffractivelv scattered light 
method are covered by Buchele [Ref. 13] and are presented 
here. 

(1) Inaccuracy of angular measurement or the limited ability 
to resolve small angles and, 

(2) Inaccuracy of the intensity measurement due to extraneous 
ligne 

57767۲1۰6 light includes all light other than scattered 
light from the particles. Some examples are scattering from 
an aperture or dirty test section windows. Refraction of the 
beam due to gas density gradients and image point broadening 
from turbulence are others. Laser speckle is also extraneous 
Ent. 


The sources of error addressed in this investigation are 


discussed in the related portions of the paper. 
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III.  EXPERIMENTAL APPARATUS 


A photograph of the apparatus is presented in +٣٥ 
A schematic is presented in Figure 2. The light scattering 
equipment was mounted on two optical benches. Comvonents for 
measurements in the exhaust plane were mounted on one bench. 
The other bench held the equipment associated with the motor 
cavity. The ligbt source was an eight (8) milliwatt Helium 
Neon laser mounted on the exhaust bench. A collimated beam 
was required so a spatial filter/collimator was used. A 
modification to.this collimator is discussed later im Ehe 
section. The collimated beam passed through a cube beam 
splitter and the second beam was diverted to a 90 degree 
prism on the other bench. The original beam continued through 
the motor exhaust plane. The other beam was routed through 
the nitrogen-purged glass windows in the motor housing. 

Each beam was then intercepted by a physical stop located 
in front of its set of receiving optics. Тһе further еше 
Stop was placed from the test section, the smaller the angle 
at which scattered light could be measured. In this 
apparatus, the stops were placed approximately 30.5 
centimeters from the exit plane of the test section. This 
allowed a minimum angle of approximately .008 radians to be 


measured. Light scattered at angles greater than ۲7 E 
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The stop served to keep the transmitted beam out of the 
measuring optics anc thus reduce extraneous light. The 
stops also improved optics alignment. This is discussed 
under calibration and evaluation of the apparatus. 

The scattered light passed through a narrow pass filter 
which admitted only light of the Helium Neon frequency. 

This filter served to reduce extraneous light from the external 
surroundings. 

An objective lens of 50 centimeter focal length was 
located behind the narrow pass filter. This lens imaged 
ао а рпосоазоае аггау the scattering profile of the "particles 
in the test section. The shadow of the beam stop was also 
imaged since the stop was between the test section and the 
Objective lens. This was a limitation which is discussed 
under calibration and evaluation. 

The photodiode arrays were the same units used by Cramer 
and Hansen [Refs. 2,12]. Each array contained 1024 silicon 
photodiodes on a single chip with 25 micron spacing. Тһе 
accompanying circuits provided a sampled and held output 
which was essentially analog except for switching transients. 
At the end of each diode scan there was a delay before the next 
scan. During this delay the diodes were reset and allowed 
to measure the intensity of the scattering peor hile aga in. 


The scanning of the diode array repeated continuously. The 
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actual sampling time of the array was about 34 miilisecome 
with a delay between scans of about 6 milliseconds. 

The 50 centimeter focul length of the objective lens 
combined with the dimensions of the diode array orovideeme 
half angle field of view of about 3 degrees for mediums of 
refractive index near one. The effective field of view was 
reduced to about 2.3 degrees for calibrations when the 
refractive index of a Plexiglas container and water was 
taken into account, 

The laser beam collimator mentioned previously at 
first produced a beam one centimeter in diameter. A lens 
in the collimator was changed to reduce the beam diameter for 
several reasons. Extraneous light would be generated if 
a large beam impinged on the aperture of the motor test 
section window. Also, if the aft beam was larger than the 
motor exhaust jet it would be refracted in the density 
(and refractive index) gradient between the exhaust and air. 

The last part of the apparatus was the rocket motor 
itself. It was the one used by Cramer and Hansen [Refs. 2,12] 


and in the present study served only for aligning the optics. 
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A. NEW CONTROLLER 

Hansen [Ref. 12] describes the major components of the 
Hewlett Packard 3054A data acquisition system. A list of the 
етиде | = relevant to this study is in Table 1. Тһе НР 85 
computer used by Hansen and Cramer was replaced with an 
HP 9836S as system controller. This newer computer has far 
more capability than the HP 85, including a choice of more 
powerful operating systems. The system used for this 
Study was Basic Extended 2.1. 

The data acquisition program written by Hansen needed 
minor modification to acquire multiple consecutive scans of 
the photodiodes. Some different I/O commands such as those 
which transfer data to the disk were also incorporated. 

The revised version of this program is listed in Appendix A. 
A general flow chart is presented in Figure 3. 

The 9836S has two internal disk drives which were used 
to store the data after acquisition. The data from both 
diode arrays was stored in the same file. The eight (8) 
scans of the motor cavity were first, followed by the four (4) 


exhaust scans. 
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B. MODIFICATIONS 

The memory capacity GE E ро вы ее 
increased so that multiple Consecutive sears of 1 
arrays could be recorded during a motor firing: THIS еи 
provide a more statistically valid measurement of particle 
size. Fluctuations o£ Scattered ао иЕе а АИ 
polydispersion need to be integrated over time or averaged 
to provide a more appropriate measurement. 

In order to fully use the memory added and make data 
management easier the data acquisition system needed to be 
modified. The memory consisted of three (3) cards, each 
with a capacity of 4096 values. The fact that this was a 
multiple of 1024 (the exact number of photodiodes) meant 
that the idle period between scans needed to be excluded 
from the data. If this was not done, one (1) less scan per 
card would have been acquired and locating the scans in 
the overall block of memory would have been more difficult. 

ТЕ was also necessary to chain two of the memory cards 
together in a way which would allow one card to be filled 
and then the other. A schematic of the data acquisition 


system is presented in Figure 4. 


Си PACING AND MEMORY CONTROL CIRCUITS 
The timing clock and blanking pulse of the photodiode 


Circuitry provided the means for pacing data acquisition. 
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(1) Memory space was fully utilized and management of the 
multiple scans made easier. 

(2) A/D conversions of the data were made exactly when a 
diode's output was on line and steady. Thus, the analog 
filter used in the previous study was no longer 
needed to suppress the switching spikes on the data line. 

Шие tO Mowing 1s a description of the signals and circuits 
used to modify the data acquisition system. All voltage 
levels were TTL. A timing diagram in Figure 5 shows the 
relations between signals. A schematic of the circuit is 

presented in Figure 6. 

The clock pulse was a positive going spike at a frequency 

Моше 30 КН. This clock Controlled all circuits of 

the photodiodes. It ran continuously, even during the blank 

period between scans when the diode output was clamped at 

zero volts. 

The blanking pulse was a signal which fell to zero at 
the beginning of each scan. It then went positive at the 
end of the scan and remained high until the next scan began. 

The clock pulse was used to drive a pulse shaper (mono- 

Stable multivibrator). This ensured that the voltage levels 

rest of the circuit would not accidentally fall‏ 50ط ت۰۰۰" 

below ЕСИ Ет 5 пеша The pulse width of the shaper 


was adjusted so that the negative going edge of each pulse 


Zo 


would occur a£ter the sw:z2cnvue transient соле саса ни 
had decayed. This negative going edge would eventually 
trigger the A/D converter to store the output of each Crees 

The blanking pulse was inverted and connected to an AND 
gate along with the pulse shaper output. The output of this 
AND gate is shown in Figure 5 as the pulse shaper signal held 
low between scans of the diode array. This was the basic 
Signal which paced data acquisition. 

This basic trigger signal was connected to an AND gate 
along with the output of the Multiprogrammers Timer Pacer 
card. In this way the trigger would not reach the A/D until 
the Timer Pacer output a pulse. This enabling pulse from 
the Timer Pacer was at least as long as the time for eight 
scans of a diode array. 

The controller programed the Timer Pacer to produce the 
pulse when the Timer Pacer received a trigger from the blanking 
pulse. Іп this way, Gata acquisition began ae the ls er oem 
a scan and no data was taken during the time between scans. 

The circuit to chain the memory cards together was 
basically an OR gate used as a negative logic AND gate. The 
end of conversion (EOC) signal of the A/D and the (FULL) Os 
of one memory card were connected to the gate. When both 
Signals went low the second card was then able to store data 


from the A/D. This arrangement is shown in the schematic of 
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the data acquisition system in Figure 4. The automatic locx- 
out feature of the memory cards when full, and the relatively 
Slow rate of data throughput made it unnecessary to control 
other handshake lines [memory card manual, Table 1]. 

The circuit was designed to handle four (4) memory cards, 
SO no modification will be necessary if one more card is 
added to the system. This would provide an additional four 


Scans of the exhaust beam. 
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V. DATA REDUCTION 


The data reduction programs written by Hansen were not 
used for this study. The new computer lended itself to 
another approach. Its memory capacity made it unnecessary to 
chain programs together and polynomial curve fitting was 
eliminated in favor of interactive graphics. Avoiding 
polynomial fits preserved the nature of the raw data so that 
one had a better feel for the parameters. The data reduction 
program "RDC" is listed in Appendix A. Figure 7 is a general 
flow chart £opyene programe 

The following is a description of the program. The user 
was first prompted for values needed to analyze a given 
data set. For example, the wavelength of the laser used and 
the index of refraction of the medium must be known for 
any data set. Next, one had the choice of reducing raw 
data scans or reviewing a reduced data file. For raw 
data one chose either the exhaust or motor cavity 
beam data. 

Raw data was plotted on the CRT and any obviously erroneous 
scans were excluded from further reduction. The valid scans 
were averaged to obtain a mean scattering profile. The mean 
intensity profile taken before particles were introduced 


was then subtracted from that taken,with particles present. 
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۰۰۰۴ or 4: Chazeocleristics O. individual photo- 
diodes and extraneous light which was independent of the 
particles. 

A symmetric moving-average-type of digital filter was 
then applied to the profile to achieve some smoothing. 

This type of digital filter was chosen for simplicity and 
because it does not have the phase lag of analog filters 
[Ref. 16]. Preserving the phase of the data was necessary 

to retain angular resolution. Another advantage of filtering 
in the software rather than hardware was that raw data 

files remained unmolested. 

The scattering profile was then analyzed using interactive 
graphics. If earlier, one chosen to review a reduced file, 
program execution began here. 

One had to normalize a scattering profile in order to 
compare it to the theoretical curves for polydispersions. 

The scattered intensity on the centerline of the beam was 
the correct value to use for normalization but was unmeasurable 
due to the beam's presence. 

The other unknown was, of course, the particle size. These 
two variables (centerline intensity for the measured profile 
and 032 for the theoretical profile of normalized intensity 
VS. (94)) were adjusted using interactive graphics until 


E. Curve tor polydiSpersions coincided with the data. Іп 


this way the mean diameter of particles was determined. 


ей 


The second reduction techn Gus Used na jp s 
application of the methoa presented by pecie ec. IM 


The equation for the polydisceme ions cung ою 
2 
L (9) = EPES И 8) 


was applied at two points of the scattering protile. шы 


gave: 


12/11 = EXP -р 7 [ (62 Ы 92) (.57 PESE 


Solving this for the diameter gave: 


dd че 2 ZU? 
0 (1 о оди / (9; 91) ] 


1 
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The computer would sweep through the data using many 
values of thetal (64) along with several angle ratios to 


determine theta2 (90 The results were presented graphically 


2s 
as particle size vs. thetal for each angle ratio (95/94). 

In actual practice the range of useable angles depends on 
the appratus, and the quality of the data. Therefore, in order 
to interpret the results one must have previously inspected 
the data. The interactive graphics routine was well suited 
to this and provided a hard сору ког ита јен 7 


After reducing a set of raw data the mean scattering 


profile was stored on disk for later review. 
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VI. CALIBRATION AND EVALUATION 


A. 70۷ ٤ ۶۹ 

The geometry of the apparatus used in the investigation 
by Cramer and Hansen is compared with that of this studv 
in Figure 8. In the previous study the transmitted beam 
was allowed to enter the receiving optics. The beam was 
focused off the diode array a few millimeters from the first 
diode. This was necessary to avoid damaging the diodes 
but introduced some uncertainty in angle measurements. 
The intense image of the beam along with scattered light from 
the receiving optics produced a high level of extraneous 
light. In the present study, stops were used to intercept 
the beam before reaching the receiving optics. These stops 
provided several advantages. A high intensity beam could 
be used while producing little extraneous light. Optics 
alignment was also improved. This reduced error in angle 
measurement. Alignment was accomplished using a neutral 
Gensity filter to reduce beam intensity and protect the 
diodes. A schematic of the apparatus 1S in Figure 2. Тһе 
laser, collimator, beam splitter and prism were positioned so 
that the beams passed through the appropriate measurement 
areas. The narrow pass filter and imaging lens were then 


ЕЕ Ее bedni entered om the centerline. Тһе‏ جج 
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photodiode array was then moved using a three-axrs nrcromc ٦ 
so that the focused beam fell on the first diode. The beam 
stops were then put in place and the neutral density rR 
removed. In this way measurements commenced exactly from 
the optical axis of the beam. 

Procedures were also refined to account for the bending 
of light rays as they passed through the walls of the particle 
container. As noted earlier, the index of refraction of the 
container and the medium containing particles affects 
scattering measurements. A Plexiglas box held the particle 
samples and a magnetic stirrer kept the samples suspended 
in water. The index of refraction of the Plexiglas and 
water combination was measured using a simple technique. 
A microscope was used to measure the ratio of actual depth 
to apparent depth for Plexiglas and water. The index of 
refraction was determined to be 1.39. This value was 
applied to the data to convert the measured scattering profile 
to that actually produced in the medium containing the 


particles. 


В. ۹٤۹۶ 

Calibration results are summarized in Table 2. Initial 
tests were done with two samples wihch were basically mono- 
dispersions of large particles. Figure 9 shows the measured 


лао Lene mses tee mec light for glass spheres ranging from 
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5750 44 microns in Giemeter. This profile was obtained by 
placing the focussed beam just far enough from the first 
diode to avoid saturation with no particles present. The 
diodes located at angles less than about .01 radians 
Saturated. The first bright ring for particles of about 40 
micron diameter was visible near .02 radians. Figure 10 
shows a profile for the same particles, illustrating use of 
the beam stops to avoid diode saturation and improve angle 
measurements. In this case, the first diode was located 
exactly on the centerline of the beam as discussed above. 

Results for a sample of 53 to 63 micron glass spheres are 
shown in Figure 11. The center lobe was nearly completely 
missed but the first two bright rings were seen near .014 
and .022 radians. Тһе first two dark rings near .01 and 
.019 radians were also seen. The method described earlier 
of setting the beam spread parameter (a8) equal to the 


zeroes and maximums of J.(a8) was used to calculate a size 


1 
Of about 58 microns. Also shown in Figure ll is the 
theoretical profile for a polvdispersion with D = 54 


D 


۲۰۰۶۹٥005۰ 

Various polydispersions of either glass spheres or 
aluminum oxide powder were then tested. These polydispersions 
consisted of fairly large particles. Results are shown in 


Figures 12 through 17. These tests showed that the apparatus 


B 


had two distinct modes ot: operation. If tne Sarticle e а 
tration was very high, or if large particles dominated the 
polydispersion, many of the diodes at Еши 8888-7 
would saturate. This left only the data at larger angles 
useable. When many diodes saturated, the theoretical curve 
given by Dobbins, et al., was used to determine size. This 
was done because this curve was valid for the larger angles 
and lower relative intensities. The curve from Buchele 
[Ref. 13] was not valid for values of the beam spread 
parameter greater than three (3). 

For low particle concentrations and/or small particles the 
data proved more accurate at the smaller angles. If no 
diodes were seen to saturate then one knew the measurement 
was in the higher intensity part of the center lobe. Here 
the curve given by Buchele was quite satisfactory for sizing. 

The smallest particles tested were five, ten, and twenty 
micron polystyrene spheres. The bright rings for these 
particles occurred at angles too large for the apparatus to 
measure. For these samples the diodes did not saturate. 
Both the Gaussian curve fit and the two angle method were 
used to obtain 032: These results were especially consistent. 
It should be noted that the two-angle method uses the equation 
for the Gaussian. ТЕ the measured profile matched the 


Gaussian exactly, Биени would be the same for any (64) 
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and angle ratio (65/84) employed. Some variations in 


52 


D due to the imperfect fit are obvious in‏ 3 0 ۶ھ 


2 
517669۲9 21 апа 23. 

A scanning electron microscope was used to photograph 
the types of particles tested. These photocraphs are shown 
uus gures 24 through 27. Equation (1) was used along with 
these photographs to calculate some of the values of 032 


in Table II. Calculations of D for the polystvrene were 


22 
arrived at using the manufacturers data on size distributions. 
The photographs generally confirm the validity of the 


technique. 


БЕ 


VII. CONCLUSIONS ۰) ۶ٰ 65 


The results of the calibration tests showed that the 
apparatus is capable of accurately measuring mean particle 
size for a broad range of mean diameters. It was found that 
the technique was most accurate if the theoretical profile 
fit or the two-angle method were applied at the smallest 
possible scattering angles. 

The rocket exhaust is likely to attenuate the beam 
somewhat, reducing the problems related to diode saturation 
at small scattering angles. Thus, measurements should be 
possible using the high intensity part of the center lobe. 
This should make data reduction less ambiguous. Actual 
testing should begin with measurements at the exhaust 
plane of the motor. These should be compared with collected 
exhaust samples to validate the use of the apparatus in 
an actual motor environment. Measurements in the motor 
cavity would then be interpreted based on the correlation 
between exhaust samples and exhaust measurements. 

It is also recommended that the index of refraction of 
the combustion gases be investigated. A literature search 
for an estimate of the index of refraction would probably 


be satisfactory. 
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TABLE I 


ELECTRONICS MANUALS 


HP Memory Cards Model 6970B Operating Manual 

HP Timer Pacer Card Model 69737A Operating Manual 
HowAsalog.to Digital Converter Card Model 69735A 
Operating Manual 

HP Users Gurde, "Using the 9826 and 9836 Computers with 
the 6942A Multiprogrammer" 

Baste Language Reference Guide with Extensions 2.0 for 


Series 200 Computers 
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1) 
20 


30 


50 
50 


80 

70 

100 
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130 
140 
150 


160 
170 
180 
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200 
210 


240 
250 
25) 
270 
280 
291 
300 
21% 
320 
338 
340 


360 
27) 
380 
390 
400 


EXC DI ACA 


PROGRAM LISTINGS 


IXXXXXXXXXX ACQDTA : FOR 98385 یعدب‎ + 
lyxxxxxxxxx ACQUIRES MULTIPLE SCANS GF BOTH DIODE ARRAYS ЫХ 
ҰЗ ХХ AND STORES THEM OX DISK عو عو لع عو وص‎ 


ы BERT HANSEN, KEN GRAHAM,KELLY HARRIS 1994 зення 
OPTION BASE 1 

COM 81$[123,02$012], T1£0201,D1$(121, D2$112] ,Address(4) ,E( 4098) 

COM A(1024) BUFFER,R(4095) ВОРРЕВ,С(4096) BUÜFFER,D(A095) RÜFFER 


OUTPUT 709; AR" IANALOG RESET 
Һ55108 88911119 TO 72310 ICLEARS THE WAKE-UP SERVICE RESUEST 
ЕКТЕК ЕН11110;041,042,002,044,045,645 IOF THE MULTIPROGRAMHMER 


MAT Address= (9) 
ENABLE INTR 7 
CLEAR 722 
98۷ ۵۷۳ 
PRINT "ENTER THE FILENAMES OF THE DATA FILES TO BE CREATED (e.g. ۸۶ 


PRINT " 12 CHARACTERS MAXIMUM, EACH" 

PRINT USING */* 

PRINT * AN EMPTY DISK MUST ВЕ IN THE LEFT DISK DRIVE ' 
PRINT '" 

PRINT ' THE DATA FILES WILL NEARLY FILL A DISK" 

PRINT USING '/* 

MASS STORAGE IS ':INTERNAL, 4,0" ! CHANGE THIS LATER IF NECESSARY 
INPUT ' INPUT FILE NAMES NOW - (FILENAME! FILENAMES) ",91%,02% 
45-1024 

Z7$=": INTERNAL , 4,1" | STRING INDICATES MASS STORAGE 


CREATE 0۵۲۱ 014, ,6 Pg 
CREATE BDAT D2$572$,0144,16 | 102 
PRINT USING #9" 


SCANS CF 1024 8 2*8 BYTES PER RECORD 
4x12/2=NUMBER OF RECORDS=5144) 


PRINT " DATA WILL BE STORED ON DISKETTE WITH FOLLOWING FILE NAMES: " 
PRINT USING "///" 

PRINT ' NO PARTICLES ---- FILENAHE = ';D1$ 

PRINT " ГАИ ЕБ جک‎ FILENAME = ",р2% 
PRINT USING '//' 

PRINT 5 15 THIS A CALIBRATION ? ENTER ’ Y ’ IF YES ' 
PRINT " ANYTHING ELSE IF NO ' 

INPUT R$ 

PRINT USING "//" 

PRINT ' BE SURE LASER IS ON" 

PRINT " PRESS (CONTINUE) WHEN READY" 

PAUSE 
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41) 
429 
430 
440 
450 
460 
470 
480 
490 
500 
519 
521 
531 
540 
55 

550 
571 
580 
591 
1 
11 
620 
630 
640 
55 
660 
679 
680 
690 
70 
718 
72 

731 
74 
750 
760 
771 
781 
791 
804 
810 
820 
834 
840 


858 
870 
880 
890 
90) 


Јин 6 l 

РА 3 + na; 

INPUT "ENTER THE THRESHOLD PRESSURE STO ІЗІСОЕК БЕ ЕВ SS 

INPUT "ENTER TIME DELAY FROM THRESHOLD PRESSURE (sec)", 18 

! TIMES ARE INTERPRETED BY THE COMPUTER IN SECONDS DOEN TO ‚001 

VO=S8/151.5 | CALIBRATION OF PRESSURE TRANSDUCER 

D3$=D1$ | NO-PARTICLES STRING NAME 

GOSUB Multiprog 

GOSUB Storedata 

IF R$="Y" THEN 1900 

| ХХ ХХХ ХХХ ХХХ УХ ХХХ ХХХ ХХХ ХХХ ХХХ ХХХ ХХХ ہو کر‎ 
Contcheck: Continuity checking 

OUTPUT 709; "AC20" ICONNECT 3456A DVM TO IGNITION FIRING LINE 

LOCAL LOCKOUT 7 IKEEPS OPR. FROM SWITCHING FRONT PANEL TO LOCAL 

7-5 | IN THOUSANDS OF OHMS 

QUTPUT 722; 'HS5M002582501L150F AR TT2QX] " IRESISTANCE MEASUREMENT 

H RESETS DVM:.SH002 SETS SERVICE REGUEST MASK WHERE 002 IS OCTAL REP. (ОР 


THE SERIAL POLL MASK EYTE;SW2 TELLS WHICH TERMINAL SWITCH IS USED;S01 
SYSTEM QUTPUT HODE ON - WAITS FOR CONTROLLER TO HANDSHAKE;L1 LOAD 
INTERNAL MEMORY ON;SQ-FUNCTION SHIFT OFF. F4-TWO WIRE CONNECTION TO DYM 
RI-AUTORAMGING; TS-SINGLE TRIGGER; Q-LOAD INTERNAL MEMORY QFF;X1-EXECUTE 
STORED PROGRAM. 

GOSUB Reading 

RS=V 

IF R8(87 THEN GOTO Ok 


Ех ххх ххх ххх хх СОНТТНИТТУ СНЕС з є ххх хх ХХХ ХХХ ХХХ 


PRINT USING "B" 

PRINT CONTINUITY CHECK ВЯ 7 

PRINT 'RECHECK BEFÜRE PROCEEDING. WHEN CHECKED, PRESS (CONTINUE J" 
BEEP 3000,.3 

WAIT „| 

ВЕЕР 100,1,0 

PAUSE 

GOTO Contcheck 


Ok:PRINT USING '8" 


DISP * СОН НОТЕ ۶۰۱۶۱ 0 ٤ 


PRINT ' BE SURE NITROGEN IS ON" 
OUTPUT ٣۳ ICONNECT DVM FOR IGNITOR VOLTAGE MEASUREMENT 
PRINT ' DUH CONNECTED TO FIRE SWITCH" 


PRINT USING '/" 
PRINT "BE SURE VISICORDER IS SET UP TO RUN ON PROPER SCALE WITH LAMP ОМ," 
PRINT USING "/" 


OUTPUT 722; "HSMO02SH2ZOSOILIFLO.O1STISOFIRIT3O" | VOLTAGE ON IGNITER 
DISP ' STANDING BY FOR IGNITION" 

PRINT ' STANDING BY FOR IGHITION' 

BEEP 2000,.1 

OUTPUT 722; "Xt" IUOLTMETER TRIGGER 
GOSUB Reading IREADS 7۶ 
RO=ABS(V) | 

IF R9(10 THEN GOTO 870 ` ! 12 VOLTS ON IGNITOR 
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910 
920 
931 
940 
759 
961 
770 
780 
798 
1000 
1010 
1028 
1030 
1040 


0-٤ Г CRD ES 
00 ٦ | CONNECT PRESSURE XDUCER 10 ۵7 
ОЏТРЏТ 722; "21" | TRIGGER VOLTRETER 

60508 ЕШ 

К9-485(У) 

IF R9(U0 THEN GOTO 950 | THRESHOLD PRESSURE 
WAIT T8 І ТІН DELAY 
O1-TIMEDATE 

IF R$O'T*. THEN 1050 ITF ACTUAL RUN THEN SKIP ОНЕ LINES 


| FOR CALIBRATIONS INTRODUCE PARTICLES AND THEN CONTINUE 
PRINT CHRS (12) 

PRINT USING '////////* 

PRINT ' INTRODUCE PARTICLES THEN PRESS КЕУ #9 ТО ТАКЕ ГАТА" 
ON KEY 9 LAREL ' ТАКЕ DATA* COTO 1060 


1030 Standby: GOTO Standby 


1060 
1070 
1080 
1090 
1180 
1110 
1120 
1130 
114) 
1150 
1150 
1170 
1180 
1190 
1280 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1500 
1310 
1328 
1330 
1340 
1350 
1360 
1370 
1380 
1390 


GOSUR Multiprog 1 TAKE PARTICLE DATA 
LOCAL 7 

D3$=D2$ IPARTICLE DATA FILE NAHE 
GOSUS Storedata 

PRINT ' LAPSED TIME: FIRE TO MULTIPROGRAHMER CALL = ":01-00 
PRINT ' ELAPSED TIME: FIRE TO MEMORY INTERRUPT 2 "02-00 
PRINT ° DATA STORED ON DISK WITH FILENAMES (";D1$;") ARD (*5D2$5;)" 
GOTO End 
Multiprog: ! 

۷ 4۶ Је, 131" 

DUTPUT 7234; "CC 577 ICLEARS THE ARM,BUSY AND ECP GF MEMORY CARDS 
POTR ۰ ۳ ICLEARS SAME РОК A TO D'S AND TIKER PACER 
Ш а 12, T ! THE (1). IS 2°S COMPLIMENT KINARY 
ИЕР ТТІ ТІГІСІ" | THE .00! 15 THE LEAST SIGNIFICANT EIT 
OUTPUT 721 "SF, УВ THEE TRE |2 To DR ИК и SIZE 


DUTPUT 723;'5F,12,3,1,.001,12,T* ! SINCE THE A TO DIS 12 ВІТ 

DUTPUT 723; "WF, 2, 1023T MF ,۸ 7 

DUTPUT 723;*F, 10, 40957 МЕ, 13, 40957" 

| SETS REFERENCE WORD FOR WHEN TO STOP TAKING DATA AND GENERATE INTERRUPT 
DUTPUT 723; "WF,2.1, 1T UF ,5. L, IT,UF, 9. ИТ, МЕ ДОД 1T? ISETS FIFO HODE 
DUTPUT 723; WF,3.1,0T,WF,3.2,0T, WF , 3.3, 0T WF, 13. 1, 0T, WF, 13.2, 0T, UE , 13. 3, 0T 


OUTPUT 723; 'WF,5.1,0T,WF 8.2, 0T WF ,6.3,0T, WF ,10. 1, 0T, WF, 10.2, 0T, 8F 10.3, 0T 


1 SETS COUNTER AND POINTERS OF 2ND HEMORY CARD IN EACH PAIR TO 0 
OUTPUT 723;AC,3T,AC,OT,AC, 10T,AC,13T°!ARHS CARDS WHICH GENERATE INTERRUPT 
OUTPUT 723; 'WF,4.2,0T,WF,4,15T* 1 TIMER PACER GIVES 1 PULSE OF 1 SEC WHEN 
| TRIGGERED BY THE BLANKING PULSE (PLENTY OF TIME FOR 8 SCANS) 
WAIT 2 
K=SPOLL(723) IWAIT FOR MEH INTERRUPT 
IF K{}64 THEN GOTO 1330 
DUTPUT 723;*WF,4.2,1T* | MAY BE UNNECESSARY TO ALTER TIMER PACER SINCE 
| MEMORY CARDS HAVE AUTOMATIC LOCKOUT BUT FOR NOW WE WILL DO IT 
O2=TIMEDATE 
ON ERROR GOTO Err trap | NEEDED TO READ ARMED CARD INTERRUPT LIST 
SEND 7,UNL MLA TALK 23 SEC 12 1 SPECIFICALLY ASKS FOR INTERRUPT LIST 
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1440 Var read: — ENTER 7;Rnodress(x) КЕЛ WHICH CA RDS INTETRUPTEZ 


1410 Memcards: PRINT ' MEMORY CARDS WHICH GENERATED IMTERRUPTS АКЕ” 
1420 PRINT ' 010754 - ';Address(x) 

1420 ЖАТ Address= (0) 

1440 OFF ERROR 

۱450000006 ,۹017 IDISSARM MEM CARDS 
1450 QUTPUT 723; "HI,2,1024T* | SET UP CARD TO BE READ 
1470 ENTER 72305;A(#) | GETS DATA FROM 1024 MEMORY BCARD 
1480 OUTPUT 723; 'HI,S 40961" 

1490 ENTER 72305; B(x) | GETS DATA FROM 4095 MEMORY BOARD 
1500 PRINT ' EXHAUST DATA ENTERED" 

15910. OUTPUT 723; ۸+ ۳ 

1500 ENTER 72105;C( x) 

1530 ۳ 

1540 ЕМЕК 72305;0(%) 

1990 ENABLE INTR 7;8 

1560 PRINT ' MOTOR DATA ENTERED" 

1570 MAT A= (-1)#A ! THE 1024 CARD IS INCLUDED RUT NOT SAVED, IT DIDN'T 
1980 MAT R- (-1)xR | PERFORM WELL, COULD ВЕ REPLACED BY A 4094 CARD, 
1592 МАТ (= (-1)*С 

1800 MAT D= (-1)*D ! DIODE VOLTAGES ARE NEGATIVE SO SIGNS ARE CHANGED 
1510 RETURN 

1520 5:1۱ 

1530  ASSIGN PDiskfile ТО 0384775 

1640 ASSIGN FBuff1 TO BUFFER C(#) ! MOTOR CAVITY 4 SCANS 
1850 ASSIGN @kuff2 TO BUFFER D(x) ! MOTOR. CAVITY 4 SCANS 
1560 ASSIGN BBuffj TO BUFFER B(x) ! EXHAUST 4 SCANS 
1570 CONTROL 8Bvff1,2;1,12768,1 ISETS RUFFER POINTERS IO FULL 
1680 CONTROL ®Buff2,3;1,32767,1 ! INTERFACE REGISTERS SECTION OF 
1590 CONTROL 8Bof f 3, 351, 327671 ! LANGUAGE MANUAL 
1700 TRANSFER @Ruffl TO @Diskfile ! ORDER OF DATA ON THE DISK IS 
1719 WAIT FOR EOT 8Diskfile ! MOTOR CAVITY--8 SCANS 
1720 TRANSFER BBoff2 TO BDiskfile ! EXHAUST -- 4 SCANS 
1730 WAIT FOR EOT @Diskfile IWAIT BECAUSE OVERLAPPING 
1740 TRANSFER @Buff3 ТО &Diskfile ITRANSFERS ARE NOT WANTED 
1750 WAIT FOR EOT @Diskfile | 

1760 ASSIGN BDiskfile TO + ICLOSE I / 0 PATHS 
1770 ASSIGN 28064 TO x 

1780 ASSIGN BRuff2 TO + 

1790 ASSIGN 850۶۶3 0 ٭‎ 

1800 RETURN 

1810 Reading :! 

1820 STATUS 7,1;A8 ICHECKING STATUS BEFORE READING 
1830 ENTER 722; ! VOLTHETER 15 A FORMALITY TO 
1840 ENABLE INTR 7;8 ICLEAR THE SERVICE REQUEST 

1850 RETURN 

1860 Err trap: IF ERRN=159 AND ERRL(Var_ read) THEN Мемсаг05 

1870 PRINT ERRM$ EVEN IF THE ERROR WAS NOT THE ONE PLANNED 
1880 ‘GOTO Hencards IFÜR PROGRAM EXECUTION CONTINUES 


1890 End: END 
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170 
180 
198 
20 


УХА АЕК ХХ ХК ХХХ КК 


عو لچ عو جع و جج ہی RDC‏ رو ود ود ڑے ےد در درد دج 
جع عو وآ PLOTS RAW DATA‏ لع جع و !ا 
XXXxXXXXXXXXXEXOXEXXX AVERAGES X XXXEXOUOHEXEX EXOEXX X‏ | 
FILTERS HERE RKA KALE KE KEK EK EX‏ رر رر ور ЗЕ‏ 


DXxxxxxiuoxxxxxxxixxx DETERMINES MEAH DIAMETER. — XX GOOGOPEDEEXOOODEGEX EXE X 
DxoaeaxxXoxoeoeG X XXX BY INTERACTIVE GRAPHICS .— еее зе × 
xxidxnxxodgoxx3xxxxxx AND THE TRO-ANGLE METHOD © #48444 KEKE eee x 
ЖЕЗ ЕЗ اصع‎ Robert Kelly Harris зе ОХ 
ЕЗУ 1984 УХ ХХХ 
VERA LK KK AKA KA AE KEELE AEE LE EEK EEK LEK ES EAE EES AEE EE EEE ENED HE EL EEE REESE 
OPTION BASE 1 

COM /Hrdaaus/ Av2(1024) 

СОМ /Gauss/ T1(1024),G(1024) ,L 

COM /Max/ M7,H5, Xt, Yt, Xn, Yn 

COM /Readata/ B,P,H,02$(1201,04$(201,2:51201,Y1(8192) BUFFER,Y2(81972) BUFFE 


COM /Tuo/ Av1(1024) M, M1,F 

DIM Scans(8),X(1024) 

COM /Plots/ P1$1201,P25(201,P35[201,P45020] 
INTEGER Graf(1:12480) RUFFER 


210 Choose: PRINT CHR$(12) 


20 
230 
Он" 
244 
250 
260 
270 
280 
290 
200 
310 
320 
330 
540 
500 
350 
370 
380 


01d-0 

PRINT ° TO LOOK AT ANY DATA FILE THE PROGRAM NEEDS SOME STARTING INFORMATI 
PRINT °° 

PRINT ° TO ACCOUNT FOR THE CHANGE IN WAVELENGTH IN THE MEDIUM" 
PRINT °° 

PRINT " ENTER THE INDEX OF REFRACTION OF THE MEDIUM" 

PRINT °° 

PRINT ' WATER = 1,33" 

PRINT ` AIR = 1,8" 

PRINT ' ESTIMATE OF EXHAUST = 1,1" 

ВЕЕР 

INPUT ' THIS VALUE ADJUSTS THE COMPARISON CURVE TO THE MEDIUM”, MI 


PRINT CHR$(CT2) 

PRINT * TO ACCOUNT FOR REFRACTION OF LIGHT AT BOUNDARIES BETWEEN" 
PRINT °° 

PRINT " THE MEDIUM AND AIR YOU MUST ENTER THE INDEX OF REFRACTION  " 
PRINT °" 
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um 
` w 
279 


PRIAT ' Cr The COMBINA:1SN Or TPz MEQIUM AhD ITS EGUNDARY' 


SUME TS 

АТОМА 7× THIS VALUE IS APPLIED DIRECTLY TO THE DATA’ 

42 PRINT "" 

430 PRINT " AIR=1,9 PLEXIGLASS k WATER = 1.39" 

440 PRINT '' 

450 PRINT "ESTIMATE QF HOTOR CAVITY COMBUSTION PRODUCTS & WINDOM - 1.22" 

450 — PRINT 55 

470 PRINT " ESTIMATE OF ROCKET АПАН 311 7٦ 

480 — INPUT M 

490 — PRINT CHR$(12) 

900 ЖИ" ENTER LASER WAVELENGTH (HeHez,5328,8rz,488)",L 

910 INPUT "ENTER THE DESIRED PLOTTING INTERVAL OF DIGDE ARRAY (2,4,6)",H 

900 ٢۵۹ им FOCUL LENGTH OF OBJECTIVE LENS 
530 — Diod-.025 Imm DIODE SPACING MAY RE .03 

940 РОВ 150 TO 1023 

jy Х(1+1)21+]1 ICREATE AN ARRAY OF DIODE NUMBERS 

250 Т1(1+1)=(1%0100)/ (ЕМ) ICREATE AN ARRAY OF THETA 

971) NEXT I 

980 PRINT CHR$(12) 

990 САП Display1(01d) 1JUST PRINTS INITIAL REMARKS ON CRT 
500 IF Qld THEN CALL Review (Avi(#)) ITHESE TWO LINES ALLOW FOR 

618 IF Qld THEN GOTO Gauss IREUIEMING REDUCED DATA FILES 

620 ВЕЕР 

530  INPUT " INPUT TWO FILENAMES NCH (NO-PART ,PART)',02$,04$ 

640 INPUT ® ENTER ‘1’ FOR MOTOR ‘2° FOR EXHAUST" ,P 

550 — Z:$z' 1INTERHAL , 4,1" ! USES LEFT DISK DRIVE FOR RAB DATA FILES 
560 CALL Readata 

670 CALL Shift 

630 KEEP 

690 Screen: |! 

700 — PRINT CHR$(12) ICLEARS SCREEN 

710  P1$-"UDLTAGE" 

72% 2$="DIODE HUMBER” | STRINGS FOR PLOTS 

730 PSS="VOLTAGE vs. DIODE" 

74) 84 P4$="RAW DATA" 

730 85-4 

7650  Xtzb4 | SET UP VALUES FOR PLOTS 

778  Ytz.04 

780 Хм=4 

790 Ya=5 

800 — CALL Plot(Y2(00) | DRAWS AXES ETC... 

810 — GSTORE Graf (4) 1 STORES GRAPHICS DISLAY JUST MADE 
820 CALL Dataplot(B,Y1(*),H) IPLOT NO-PARTICLES DATA 

830 РЕТ! NO-PARTICLES PLOTTED" 

840 PRINT "" 

850 PRINT ' HOW MANY SCANS SEEM TO EE GOOD 7?" 
860 INPUT J 
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373 
35% 
MMA? 
890 
700 
910 
920 
930 
740 
950 
750 
970 
984 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1050 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 


PAINT °° 

PRINT ' WHICH SCANS ARE GOOD 7, ЕАО 0 
НАСТ СОЯНА 15 REQUIRED. A YOU HAVE TO HIT CONT.. TWICE 

INPUT Scans{*) 

PRINT CHRS( 12) 

PRINT USING “////////7" 

PRINT ' AVERAGING THE SELECTED SCANS" 

GRAPHICS OFF 

CALL Average(J,Scans(*) ,Y1(*) АУ1(%)) 

PRINT USING '8' 

GLOAD Grafix) !LOADS GRAPHICS ARRAY RATHER THAN WASTE TIKE RE-DRAWING 

GRAPHICS ON 

CALL Resul t(Avi(*) ,X(*),H) ! PLOTS 1024 ELEMENT ARRAYS 

PRINT USING '//* 

PRINT ° Average Intensity No-Particles' 

QN KEY 0 LABEL 'AVERAGET" GOTO Screen 

ON KEY 1 LABEL "PLOT-PARTICLES" GOTO 1059 

PRINT ' PRESS KEY $ 0 TO RE-AVERACE OR 3 1 TO CONTINUE” 


Standby: | MANUAL CALLS THIS INTERRUPT DRIVEN PROGRAMMING 


GOTO Standoy ! LOOPS, WAITING FOR USER TO DECIDE 
PRINT USING "В" 
OFF KEY 1 IHELPS AVOID CONFUSION BY CLEARING THAT BOX 
CALL Dataplot(B,Y2(x),H) PLOT PARTICLES DATA 
PRINT USING "/" 
PRINT ' PARTICLE DATA PLOTTED" 
PRINT * FOR A HARD COPY GF THIS RAW DATA PRESS KEY # 6” 
PRINT °° 
PRINT " OR TO CONTINUE PRESS KEY # 1° 
ON KEY & LABEL "HARD & RAW" COTO Raw 
ON KEY 1 LABEL "CONTINUE * GOTO Select 
GOTO Standby 


1170 Raw: CALL. Plot( 200,1) | THE ONE (1) IS AN OPTIONAL ۶۲ 


1180 


CALL Dataplot(B,Y2(9),H) | WHICH IS USED TO GET HARD COPIES 


1190 Select:! 


«1280 
1210 
122 

1230 
1240 
1250 
1250 


OFF KEY 6 
OFF KEY 1 
PRINT CHR$12) 
PRINT USING *////" 
PRINT ' HOW MANY SCANS SEEM TO ВЕ GOOD ?' 
INPUT J 
PRINT ° WHICH SCANS ARE GOOD ?...ie...1,2,4,5,7, INCLUDE 


LAST СОННА" 


127) 
1280 
1290 
1300 
1310 
1320 
1330 
1340 


INPUT Scans{*) 

PRINT CH8R$(12) 

PRINT USING 0074۰ 

GRAPHICS OFF 

PRIN S AVERAGING THE SELECTED SCANS" 
CALL AveragelJ,Scans() , Y2(x) ,Av2(3)) 

PRINT USING "8" 

PASZ'AVERAGED SCANS” 


al 


1350 
| 350 
1370 
1380 
1390 
1400 
1410 
1420 
1430 


ПАН аи ТІ 
GRAPHICS ON 
ТОЕ бга? (ж) 


CALL Result(Av2(%) ,X(%),H) 


PRINT USING "//* 
PRINT ' 


Average Scattered Intensity” 


ON KEY 2 LABEL "AVERAGE2" COTO 1080 
ON KEY 3 LABEL "SUBTRACT" GOTO 1450 


PRINT " 


CONTINUE” 


1440 
1450 
1450 
1470 
1480 
1490 
1500 
1510 


GOTO Standby 

PRINT USING "8" 

DFF KEY 3 

MAT Avi= Av2-Avi 

GLOAD Graf(x) 

CALL Result(Avi(t X), XC), H) 
PRINT USING. "///* 

PRINT ' 


0-1۳ 


1520 
1510 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1520 
1530 
164) 
1550 
1550 
1570 
1680 
1690 
1200 
1718 
172) 
1730 


PRINT ° 

PRINT ' 

QN KEY b LABEL ° 
OH KEY 1 LABEL " 
GOTO Standby 
CALL Plot(Avi(®), 1) 

CALL ResulttAvi G0, XCO,H) 
PLOTTER IS 3, "INTERNAL" 
PRINT USING "//" 

PRINT ' 

PRINT " 

PRINT " 

PRINT " 

INPUT " 

IF Fil=0 THEN Gauss 
P4$="FILTERED DATA” 

CALL Filter (Av1(#),Fil) 
28 

CALL Plot(Avitx)) 

IF Old THEN HAT X= ТІ 

CALL Result(Avl(*) ,X(*#),H) 
PRINT °" 


FILTER 


0 ۹٦ 


1740 
TER" 
1750 
1760 
1770 
1780 
1790 
1800 
1810 


PRINT " 


PRINT ° 
PRINT ° 
PRINT ° 
PRINT " 
PRINT ° 


HIT KEY # 2 TO RE-AVERAGE OR KEY + 5 10 


ISURTRACTS NO-PARTICLES FROM PARTICLES 


Plot of the Difference Between Particles and N 


KEY # & FOR HARD СОРУ" 
KEY #1 TO CONTINUE " 


HARD AVERAGE” GOTO 1570 


" GOTO 1400 


ENTER THE NUMBER OF TIMES YOU WISH TO APPLY” 
THE DIGITAL FILTER FOR SMOOTHING" 
EXAMPLE RRERKARX 10" 
TAKES ABOUT 1,5 MINUTES" 
ZERO 15 ALSO ACCEPTABLE” ,Fil 


Plot of the Difference Between Particles and N 
AFTER APPLICATION OF A DIGITAL FIL 


for a HARD COPY ® 
PREPARE the PLOTTER and” 
PRESS KEY # 8” 
08" 
PRESS KEY # 9 TO NORMALIZE" 


OX KEY 8 LABEL "HARD FILTERED" GOTO 1839 


ON KEY 9 LABEL * 


NORMALIZE 


" 6010 Gauss 
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192) 
1830 
1840 
1850 
1850 
1870 
1880 
1890 
1900 
1910 
1929 
1930 
1940 
1950 
1750 
1970 
1980 
1990 
2000 
2010 
2320 
2030 
2040 
2050 
2050 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2159 
2170 
2180 
2190 
2208 


224) 


r 
r 


го го го го га га го 
tJ od PO ГӘ гә ҒӘ ғо 


- => ~O CO чу се cn 
=> = <> c => c < 


GOTO Standby 

8١۱ 4 04047 |! The | is far hard copy 
CALL Resuvlt(Avl(e) (ж) И) 

FRETS 


Gauss: OFF KEY 2 IKEY 1 OPTION IS STILL VALID 


OFF KEY 8 
ON KEY 1 LABEL ° FILTER * GOTO 1600 ! ALLOWS ONE TO SMOOTH OLD DATA 
| WHEN REVIEWING REDUCED FILES 
PRINT CHR$(12) 
lAv2 Array is Normalized BUT Av1(*) 15 Still Saved For Re-werk if Needed 
MAT Ave= Avi — ! TO REGIN WITH THE ARRAY IS ASSUMED TO BE NORMALIZED 
! AND TRE USER CHANGES THIS WITH INTERACTIVE GRAPHICS 
! IN THE SURROUTINE "Compare 

P1$="NORMALIZED INTENSITY" 
P2$="THETA (rad)" ISTRINGS FOR PLOTS 
ЕНЕЛ 5 HE roe 
PES CURVE FIT RESULTS" 
805 ISET UP VALUES 
Xt=.004 IFOR PLOTS 
Yt=.1 

Хи=2 
Үнгг 
CALL Plot(Ave(+#)) 
CALL Result(Av2(%) ,T1(4) ,H) IPLOT OF NORKALIZED INTENSITY vs. THETA 
СТОКЕ Graf (¥) | SAVES THE SCREEN IMAGE 
CALL Cospare(Avi G0 ,H, M, H1, Graf GO) 

GRAPHICS OFF 

ON KEY 4 LABEL "OTHER ARRAY" COTO Choose 

ON KEY 5 LABEL °" TWO-ANGLE" GOTO Buchele 

ОН KEY 9 LABEL * NORMALIZE" GOTO Gauss 

ON KEY & LABEL "STORE DATA" GOTO Storedata 

ON KEY 7 LABEL " QUIT ° GOTO Quit 

CALL Menuel 1! PRINTS OPTIONS ON THE SCREEN 
GOTO Standby 


Buchele: |! 


CALL Twoangle 
CALL Menvel 
GOTO Standby 


Storedata:! 


CALL Store 
CALL Menvel 
GOTO Standby 


Quit: | 


END 
SUB Average(d ,Scans(*),Y(*) ,Ay(*)) AVERAGES SELECTED SCANS 
MAT Av= (0) ! INITIALIZES THE ARRAY LOCAL 19 THIS ROUTINE 
РОК 151 T0 J STEP 1 
K-(Scans(D-1)x102441 ` | THIS COUNTER IS THE REGINNING 
L-Scans( 211024 ! AND THIS ONE THE END OF BLOCKS 
ҒОК П-К TO L ! OF 1924 IN THE OVERALL DATA 
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Цена ۰۱ ! ۱۴ )  + +9 oc s 


S = 

2330 ВУ) ЕУ СПОЈЕНЕ | | AND 1024 

2340 NEXT Ti 

2530 NEAT Í 

2360 MAT Av= Av/(J) HOeexTHE. AVERAGE (9G GG! 

2370 SUREND 

2380 SUP Dataplot(B,YOO А) !8 IS 4 OR B(THE NUMBER СР SCANS) 


2390 COM /Max/ 87,5, ۷۸ 
2480 LDIR 0 


2418 LORG 4 

2420 FOR J=1 10 8 IEACH SCAN 

2430 МОЧЕ 50047x50,47- . 03 | SEE NEXT LINE 

044$ LABEL J; | HELPS KEEP TRACK QF EACH SCAN AS IT APPEARS 
2450 LINE TYPE 1 

2460 K=(J-1)#1024+1 IBEGINNING OF EACH SCAM AND 

2470 L=J*1024 ITHE END WITHIN THE TOTAL BLOCK 
2480 MOVE 1, YOO |! MOVE TO THE FIRST POINT 

2490 FOR I-K TO L STEP H 

2300 X=I-K+] | THIS GIVES 1 TO 1224 FOR ABSCISSA 
2510 PLOT X (I) 

2520 NEXT I 

0530 НЕХТ ) 

2940 PENUP 

2550 1 

255 SUB Result(Iav(,DGO,H) ITHE AVERAGE INTERSITY IS PLOTTED 
2579 CLIP ON 

2580 MOVE D(t), ач (1) 

25970 FOR I=1 10 1024 STEP H 

2600 ОКАВ 001), Таи (1) 

2610 NEXT I 

2620 PENUP 

263 SUBEND 

2540 SUB Plot(Y¥(*) OPTIONAL Hard) !IF THE OPTIONAL(Hard)IS RECEIVED 
2550 COM /Max/ M7,MS, Xt, Yt, Xn, Yn ITHE FIGURE GOES TO THE PLOTTER 
2660 COM /Plots/ P1$020),P23$(20],P 3$(20] ,P4$(20] 

2670 = M7=HAX(Y(#) )#1,1 ! A SCALING VARTARLE 

0680 16۴,1 | ANOTHER SCALING VARIABLE 

2690 PRINT °° 

0700  СІНІТ IGINIT IS JUST GCOD PRACTICE SO YOU KNOW 
2718 IWHERE YOU ARE BEGINNING 

2720 SELECT NPAR ITHIS DETECTS IF THE HARD COPY IS DESIRED 
2730 CASE 1 IIT COULD FE DONE WITH IF THEN LOGIC 

2740 PLOTTER IS 3, "INTERNAL" ВЫТ IS PRESENTED FOR FAMILIARIZATION 
2730 CASE 2 ISINCE IT IS MORE POWERFUL IN COMPLEX 
2760 PLOTTER IS 705,'HPGL" ۷۵ 

0770 END SELECT 

2780 GRAPHICS ON ! TO BE APLE TO SEE THE PLOT 
2778 6۹۵ | 

2800 DES 

2810 LDIR 0 
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20 
2% 
4) 
59 
60 
70 
880 
89) 
200 
710 
920 
759 
940 
950 
2961 
297$ 
2980 
299) 
3000 
3010 
3020 
3030 
3040 
3050 
3050 
3070 
3088 
3090 
3103 
3118 
3120 
3130 
3140 
3150 
3160 
3170 
4180 
3190 


ә Го ГӘ ГО ГӘ ГО ГӘ го ГО Го ГО гз ГО ГО ГО 


3200 
3210 
3008 
222 
324) 
3250 
3250 
3270 
3280 
3290 


LORS 5 
HOVE 75,95 
LABEL USING "K';P4$ 
MOVE 75,90 
LABEL P3$ 
MOVE 75,2 
LABEL P2$ 
LDIR 90 
MOVE 18,50 
LABEL P1$ 
VIEWPORT 30,125,30,85 
FRAME 
WINDOW 0,M5,-M5, M7 
AXES Xt ,H5,0,0, X», 1,5 
LDIR $ 
CSIZE 2,.5 
LORE 8 
CLIP OFF 
FÜR Iz-H& TO M7 STEP 87/10 
MOVE 0,1 
LABEL USING "4,HD,DDDD*s1 
NEXT I 
CSIZE 3,.6 
LDIR 90 
LORE 8 
FOR Iz0 TO M5 STEP Xt 
МШЕ 1,-М5 
LABEL USING "i,K"I 
NEXT I 
PENUP 
SUBEND 


SUR Conparetàv1 OG, H, H, AT, INTEGER. Graf G0) 


COK /Gauss/ T1(*),G(*),L 
COM /Hrdgaus/ Av2(X) 
COM /Max/ M?,H5,Xt , Yt Xn. T^ 


| SCREEN UNITS FOR MARGINS 
I DRAWS A BOX 
|I EXTENT OF X AND Y 


INUMBER THE Y AXIS 


INUMRER THE X AXIS 


DIM T(1:1024),51(1:17),D0bbins(1:17) ,T1d (1:17) 

ШЕН 4055355559 6, 6.59,7.,7.9,8, , 9.3,2.,9, 9, 18. lea, tl, 
DATA „0596, 039, 0206, „014, 0106, 0081, .00605,.005, 00374, 0032, 00248, 
.0022,.00185,.0016,.00115,.0012,.001 


READ At(X) 

READ Dobbins(xX) 
OFF KEY 

ON KEY 3 LABEL 'MENUE" GOTO Subexit 
PLOTTER IS 3,"INTERNAL' 
Нагд= 

Септег Цпе-| 

0:08, 

CSIZE 4,.6 

PRINT USING ОДА 


Change: 


5 


IGETS RID OF ALL LABELS ON KEYS 


| IN CASE A HARD COPY WAS JUST MADE 

! (0) SO ONE DOESN'T EXIT TOO SOON 

| THE INITIAL NORMALIZING VALUE 

| INITIAL PATICLE MEAN DIAMETER IN MICRONS 


3300 PRINT ' STAND BY FUR CURVES 
3310 МОЧЕ 0,1 

342) CLIP ON 

3330 ConzPIsDxM1/L ГА CONSTANT 

3340 FOR I=1 10 1024 STEP H 

335) 01 0 08 ITheta-bar for some Фбаг(лоте 81) 
3350 6) ЕР ои Gaussian curve for this Dbar 
1 IF T(D23,0 THEN 3400 IGAUSSTAN NO COUD FOR Theta-bar)J 
4380 DRAW T1(D,60D IGayssian curve vs, theta 

3399 NEXT I 

4400 MAT Tid= At/(Con) (THETA FOR POINTS ON DOBBINS’ CUR 
VE 

3410 FOR I=! TO 17 

4420 DRAW Tid (I), Dobbins(I) 

4430 NEXT I 

3440 PENUP 

4450 IF Hard THEN SUBEXIT IEXITS ROUTINE IF A HARD COPY WAS JUST MADE 
3460 PRINT USING "/////////* 

4420 PRINT " USE THE KNOB TO VARY THE PARTICLE 
Sees 

3480 PRINT " OR HIT KEY # 6 FOR HARD COPY" 
4490 PRINT ° * 

2500 PRINT ° KEY $ 9 ALLONS YOU ТО" 
3510 PRINT " RE-NORMALIZE" 

352% PRINT °" ° 

3530 PRINT " HIT KEY $ 3910 бер ОЛ” 
3940 PRINT USING '//////* 

3550 QN KEY & LABEL 'HARD COPY* GOTO Hardgauss 

3560 ON KEY 9 LABEL “NORMALIZE” GOTO Divide 

3970 ON KNOB .05 GOTO Pulse 1(,05) IS TIME INTERVAL IN WHICH 
J990 Wait: GOTO Wait ! PULSES FROM THE KNOB ARE 

3390 Pulse: PRINT 05146 "9" ! COUNTED AND THIS NUMBER IS 
3600 Strngl$z'Mean Diameter ' | USED BY THE INTERRACTIVE 

3610 trna2$-' Microns' | GRAPHICS TO VARY THE PARTICLE 
352% Coun t=KNORX ! SIZE AND PLOT THE ASSOCIATED 
3530 D-DROUND(D4Count/15,2) ! GAUSSIAN APPROXIMATIONS QF 
3540 CLOAD Graf(x) | SCATTERING PROFILES 

3550 LDIR 8 

2550 LORG 8 

3670 MOVE AS, .B*47 

3680 LABEL Strngi$&VAL$(D)aStrng2s 

3590 GOTO Change 

3700 ٦٥۴٥۰۱ 

3710 PRINT ° PREPARE the PLOTTER" 
3720 PRINT " PRESS CONTINUE fer" 
3730 PRINT °" HARD COPY" 

3740 PAUSE 

3739 Нагд=] ! 50 THAT SUBEXIT OCCURS AFTER HARD COPY 
4750 CALL P1ot(Av2(3),1) | pexHARD COPY 

3770 CALL Result(Av2(*) ,T1(#) ,H) 
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5780 
4790 
3808 
3810 
3820 
3830 
3840 
3850 
3850 
3870 
3880 
3890 
3900 
3910 
3920 
3930 
3940 
2950 
3960 
3970 
3980 
3990 
4000 
4010 
4020 
4038 
4049 
4050 
4050 
4070 
4080 
4090 
4100 
4118 
4120 
4138 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4261 
4270 


PRINT USIO ipi 
LORG 9 
MOVE KS, .S%H7 
LDIR 0 
۸2٤ح‎ 5 6۹ 
MOVE 0,1 
GOTO Change 
Divide:! 
PEN -1 IERASES LAST PLOT 
MOVE 0,0 
FOR I-1 TO 1024 STEP H 
DRAW T1(1) ,Av2(]) 
NEXT I 
PRINT CHR$(12) 
PRINT USING '//////* 
PRINT ° USE THE KNOR TO RE-NORMALIZE" 
PRINT ' KEY # 4 WILL RE-PLOT USING" 
PRINT ° UPDATED SCALE" 
ON KEY 4 LABEL " RE-PLOT” COTO Replot 
ON KNOB „05 GOTO Vary 
COTO Wait 
Vary:PRINT CHR$(CI2) 
EER !ERASES LAST FLOT 
MONE 0,0 | 
FÜR 121 TO 1024 STEP H 
DRAM T1 CI), Av2t1) 
NEXT I 
PEN 1 
C2-KNGBX 
Centerline-Centerline*C2/120 IKNOB USED TO VARY NORMALIZING VALUE 
MAT Av2= Avi/(Centerline) INEW NORMALIZED DATA ARRAY 
GOTO 3870 IPLOTS THE NEW NORMALIZED ARRAY 
Replot: | 
PRINT CHR$(12) 
CALL Plet(Av2(X)) 
CALL Result(Av20) , T1303), H) 
GSTORE Graf(X) 
GOTO Change 
51۱ | 
SUBEND 
SUB Menvel 
PRINT USING ۰۳ 
PRINT ° YOU CAN RE-NORMALIZE' 
PRINT وپ بب‎ EK EEK RE KE KI پر روغ یرجہ )رر‎ XXIX XXX ЕЕЕ EC X X" 
PRINT ° YOU CAN RE-AVERAGE (Mew Data Only)" 
PRINT وع وع ہو وع عو وپ رو ور وع ہپ‎ КИК ХЕХЕ ЖЕ ۴ع یع‎ 
PRINT ' PRESS KEY # 4 TO LOOK AT OTHER DATA‘ 
PRINT ° OLD OR HEY ° 
PRINT °? MOTOR / EXHAUST” 
PRINT XXX HK HEH LH EH EK KK EEK EEK LEK EE EERE EER RK HER EKER X 
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4280 
4290 
4300 
4310 
4320 
4250 
4340 
4350 
4350 
4370 
4380 
4390 
4400 
4410 
4420 
4430 


PRINT ' PRESS Ke r lDN 
PRINDS FOR THE ’ TWO-ANGLE ’ METHOD" 
PRINT * OF PARTICLE SIZING” 
PRINT "'XYXXYXYXXYXXXYYXYYYKXXXKYYAYYYKYYYKXXYYXXXYYXYYEXXYEYYXXYYKXKKXX? 
PRINT ° TO STORE THE REDUCED DATA PRESS KEY # & е 
РАТНТ "ХУЖЖХХУХУУУ ЖУК ХХХ УХ УХУ ХЕ ХЕ ХХХ * 
PRINT ' TO END THIS SESSION PRESS KEY # 7” 
SUBEND 
SUR Displayi (01d) 
PRINT ° TO REDUCE NEY DATA PRESS KEY $ 1? 
PRINT '' 
PRINT ° TO REVIEW PREVIOUSLY REDUCED DATA PRESS KEY 3 2' 
PRINT" ° 


ОН KEY 1 LABEL ° NEW DATA" GOTO New 
ON KEY 2 LABEL " OLD DATA" GOTO Review 


Wait: GOTO Wait 


4449 New: PRINT CHRS(12) 


4450 
4459 
4470 
4480 
4490 
4500 
4510 
4520 
4538 
454} 


PRINT ° PUT THE DISK IN THE LEFT DRIVE IF IT IS NOT ALREADY" 
0 ™' 

PRINT ° ENTER THE NAMES OF THE FIRST AND SECOND FILES.” 
PRINT ** 

PRINT * EACH FILE HAS DATA FROM BOTH DIODE ARRAYS,” 

PRINT "'' 

PRINT * ` FIRST IS NO-PARTICLES & NEXT IS PARTICLES ----D1$,D25" 
PRINT "°" 

0180 

SUBEX IT 


4538 Review: 014=1 THIS VARIABLE IS PASSED TO THE MAIN PROGRAM TO INDICATE 


4560 
457% 
4580 
4590 
4600 
4619 
4628 
4530 
4540 
4550 
4558 
4570 
4080 
4590 
4700 
4710 
4708 
4730 
4740 
4750 
4750 
4770 


ІТНАТ ТНЕ DATA ТО RE READ IS ONE (1) SCAN AND THAT 
ITHE AVERAGING ROUTINES WILL NOT RE USED 
SUBEND 
SUR Readata 
СОМ /Readata/ B,P,H,03$0201,04$1201,22$020] , Y1(4) RUFFER,Y2(x) BUFFER 


TUIS ! FILE POINTER VARIARLE 
B=8 ! NUMBER OF SCANS IN MOTOR DATA 
R 1-559385 'NUMBER OF BYTES OF MOTOR DATA 
IF P=2 THEN Xyz=4097 'RECORD $ WHERE EXHAUST DATA BEGINS 
IF P=2 THEN R1=32768 INUMBER OF BYTES GF EXHAUST DATA 
IF P=2 THEN B=4 'NUMBER OF SCANS IN EXHAUST DATA 
PRINT °" 
PRINT ° READING DATA FROM FILE ON DISK" 

ASSIGN Disk! TO 0344775 

ASSIGN 8Disk2 TO 04$42:$ 0 0 ٤ 


ASSIGN ۵8۶۴1 ТО 80۶۶٢۶ Y1(%) 
ASSIGN 8Buff2 TO BUFFER Y2(#) 
CONTROL $Diskl,5:Xyz 


TRANSFER &Diskl TO PBuff1;COUNT R1 IREADS NO-PARTICLE DATA 
WAIT FOR ЕОТ 801511 ! NOT AN OVERLAPPING TRANSFER 
CONTROL @Disk2,5;Xyz Xyz SETS DISK FILE POINTER TO 


{EITHER MOTOR OR EXHAUST DATA 
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478) 
479) 
4990 
4810 
4820 
4830 
4840 
485 
4850 
4870 
4880 
4370 
4900 
4910 
4920 
4930 
4940 
4750 
4950 
4970 
4980 
4990 
3000 
3010 
9028 
5030 
3040 
9050 
2050 
5070 
5080 
2090 
2100 
2110 
2120 
5140 
5140 
3150 
2150 
9170 
9180 
0190 
5000 
0010 
5220 
5230 
9240 
2250 
9060 
2270 


TRANSFER SDisk2 TO #buff2 COUNT RI 


WAIT FOR EQT SDiske 
ASSIGN @Diski TO * 
ASSIGN @&Disk? 10 x 
ASSIGN @Ruffl TO x 
ASSIGN @Buff2 TO x 


SUBEND 
SUB P1ot2(E,D GO ,Х) 


DIM 11 1 01 
CINIT 
IF X=1 THEN PLOTTER IS 703, “HPCL" 
DEG 
GRAPHICS ON 
VIEWPORT 35, 125,30,85 
Max=10%INT( (MAX(D(#))410)/10) 
WINDOWS D(1,2), D(E, 22,0, Max 
б-(%(Е,2)-0(1,2))/(Е-19%4 
F=(INT(E/4)-1) #4 
IF F=0 THEN F=4 
CLIP D(WA,2)-G, DCF ,2) 256,0, Max 
CRID G,5,D(4,2)-C, 8 
CEIR OFF 
LORC 8 
LDIR 90 
CSIEE 90 
FOR 1=4-T0 E STEP 4 

MOVE D(I,2),0 

LABEL USING '.DDDD';D(I,2) 
NEXT I 
LDIR 0 
LORC 8 
FOR J=10 TO Max STEP 10 

MOVE D(4,2)-6,1 

LABEL USING '$,K'*;I 
NEXT I 
PSI7E 6) 5 
Titles="TWO-ANGLE METHOD" 
Titlel$-'For Various Angle Ratios" 
oize$z'SIE (microns)' 
Thta$="THETA  (rad)" 
бури! |" 
LDIR 90 
LORC 5 
B=D(4 ,2)-G-(D(E,2)-D(1,2))/18 
MOVE B,Max/2 
САВЕ. 51206 
LDIR 8 
Az(DCE,2)4D(1,2))/2 
HOVE A,-Hax/4 
LABEL Thta$ . 


po) 


I READS PARTICLE DATA 


JUST GOOD PRACTICE TO CLOSE 
11/0 РАТН5 


IDEGREES fer LABEL DIRECTION 


!AN X GRID LINE EVERY 4th POINT 


MAKES GRID UNIFORM 


IPUTS NUMBERS ON X AXIS 


NUMBER Y AXIS 


! STRINGS FOR PLOTS 


5290 LORS 3 

5290 ٦0۷۴ ۵--۵۹۵ 
3300 LABEL Sub$ 

531) LORG 5 

5300 MOVE A ۹٤ص‎ 
5330 LABEL Titles 
340 CSIZE 4.5 

3390 MOVE А, Мах! |, 05 
5350 LABEL Titlel$ 
2370 PENUP 

9389 SUBEND 

9390 SUB Distribution(D(*) , Tratio,£) 
3400 LINE TYPE | 
3410 CLIP QN 


5400 Tratio$zVALS(Tratio) ! ANGLE RATIO 
9439 HOUE DCT, 2), DCT, DD 

9440 FOR I=1 TOE ITHIS SUBROUTINE PLOTS THE 
9450 DRAM DCI,2) , DCI, 1) IPARTICLE SIZES DERIVED USING 
9450 НЕХТ І IVARTOUS ANGLE RATIOS APPLIED 
5470 LORG 4 ITO THE DATA OVER A RANGE OF 
3480 BL E 1 ! ANGLES 


5490 © CSIZE 4,.3 

1980 ПОЕ ЕЕ ٥۶ 

5510 LABEL Tratio$ 

5920 PENUP 

9548 1 

33940 SUB Tuoangle 

7990  DO900000000000090C0G0DQCQOOGOOOGQQQOOOQOGOQQQ€ECEQOQXOQOQOCEEXOODOOOGUEEXEXXGEXEE 


9560 | ІГзевевевеевеено а SURPROGRAM 'TWOANGLE' 221223532225 532255%- 
5570  UeexxeEpeee ect PARTICLE SIZING BY XXX XXX XX KK XX 
5580  beeeeeceoeeExxx MEASURING THE RATIO XXXERXXOEOEKXXOOGOEX XXX 
g590  heeeeeeeoexo xx OF INTENSITY AT ЖХХЖХХ ЖЖЖ ЖЖЖ ЖЖЖЖ %% 
5600 | Рзвевевеєне ою в THO ANGLES ХХ ХХХ 


SALO LEERE E E E E E EE E IEE E E IE IEEE E E E EE E E E EE EEEE EE E A K EEE EE E E EE EE EEE E E EE 
J620 OPTION ВАЗЕ 1 
5630 COM /Two/ Avi(1024) ММ1,Е 


9640 СОМ /Gauss/ T1(1024),6(1024),L ITHE GAUSS JS NOT USED HERE 
5650 IBUT THE COM BLOCK HAS THETA 
5650 DIM D(200,2 LAND WAVELENGTH 


J670 PRINT CHR$(12) 

9680 PRINT USING "///////" 

5690 PRINT ' THIS SUBPROGRAM USES THE TWOÜ-ANCLE METHOD DESCRIBED BY BUCHELE' 
4/08 PRINT ' " 
4710 PRINT 


р TO CALCULATE PATICLE SIZE FOR VARIOUS ANGLE RATIOS AND ANGLES." 
4/20 PRINT ' " 


3730 PRINT IT IS HOPED THAT THE CURVES WHICH RESULT WILL SHOW " 
3740 PRINT " ' 
W700 PRINT " WHICH SIZE IS THE MOST PROBABLE." 


37560 PRINT " ° 
3770 PRINT "AFTER NOTING FROM THE RAW DATA WHICH ANGLES CONTAIN THE CENTER LORE" 
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۵ ۶8۷۷ھ 

3/90 PRINT ° ВЕ мае рю ГЕ ОМА ЕТА саан Е" 

“$00 PRINT " E SMA LLEST ANGLE RATIO, ARD” 

3810 PRINT ' THE STEP SIZE BETWEEN ANGLE RATIOS" 

9820 PRINT " YOU WISH TO EXPLORE" 

3830 PRINT ' ۸۵۳۶۲۲۱۴۴۴۴۴۴ „Ме ен“ 

3840 INPUT ' ENTER thetal,theta-ratio,theta-step', Q,A, B 

2850 DER KEY 

3860 ا‎ IGRAPH WILL APPEAR ON SCREEN 
3870 Begin: IRUNNING CONTINUES HERE WHEN A HARD COPY IS DESIRED 
2880 IVan De Hulst and Gumprech & Sleepevich explain that the chanae in 
5870 lwavelength of the beam 15 accounted for by dividing by the index 
3200 lot refraction of the medium, #txxžxxxx L =L/M1 эн» 

3910 Cz(L/4317.57/PT)^2 حم ا‎ -see page 19 ef nasa tech paper 2156 
2920 РОК 121 TO 1024 Ito see this is a convenient constant 
5930 IF GCT1080. THEN. $950 IFINDS POSITION OF MINIMUN ANLGE 
3940 NEXT N 

9990 FOR TratiozA TO 3 STEP B IVARIDQUS ANGLE RATIOS 

9950 FOR J=X ТО 1010/Тгатіо STEP 10 ISETS THE RANGE OF POSSIBLE 
3570 | ANGLES (THETA 1) 

3980 Та = И (Ј) 

3990 Thta2zTI (JaTr atio) ! THETA 2 

5000 DeltathetazThta2*2-Thtaí^2 

6018 115] 

5020 I2-0 

5030 РОК 1=Ј-5 ТО Ј+5 ITHESE 2 LOQPS AVERAGE A FEW 
5040 Пере (1) !INTENSITIES IN THE HÜPE 
5050 НЕХТ 1 IGF A STEADIER CALCULATION 
6060 11=11/11 

5070 РОК 12191 (Је Га! 10-5) TO INT(J*Tratiots) 

6080 12=12+4%1(1) | DONE HERE FOR THETA? FOR 
5090 МЕХТ 1 ! THE GIVEN ANGLE RATIO 

6100 12=]2/11 

6110 E=(J-N)/10+1 ITHIS IS A COUNTER FOR THE ARRAY CONTAINING 
6120 IPARTICLE SIZE AND THETA] FOR THE GIVEN 
6139 LANGLE RATIO 

6140 IF I1(12 OR 1120 OR I2(20 THEN D(E,12z0 !ALLGUAMCE FOR IF THE ФАТА 
6150 IF I1(I2 O8 I1z0 QR I2(20 THEN GOTO Xcomp 115 NOT WELL BEHAVED 
6160 D(E,1)2S5GR(-C/DeltathetasxLOGCI2/I1)) ! ۲8۷١۵ ANGLE METHOD 
6170 IVALUE OF DIAMETER based on INTENSITY RATIO 
6188 For a given ANGLE RATIO and ANGLE THETA1 
6190 Хсомр: D(E,2)=Thtal 

6200 Spar-zPI*D(E,1)*M1/L ITHIS IS PI¥D/LAMBDA -- THE SIZE PARAMETER 
5210 83آ‎ ۵ ITHETA BAR FOR THE LARGER ANGLE 
5220 IF Tbar?i THEN J-1010/Tratio ITHIS ENDS THE DO LOGP FOR THIS 
6230 IF E21 AND Tbar)3 THEN 5330 !ANGLE RATIO SINCE THE GAUSSIAN 
6240 the above line ends all calculation if 115 NOT VALID WHEN Tber ) 3 
6250 ‘the first element (E=1)failed the test 

6250 NEXT J | 

5270 | FIRST TIHE THROUGH-- D HAS THE MOST ELEMENTS IT WILL HAVE 
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5296 
5290 
6399 
5310 
6320 
6330 
5340 
6350 
6360 
6370 
6380 
6390 
5400 
5410 
6420 
6430 
6440 
5450 
5460 
6470 
5480 
5490 
6500 
5510 
5520 
6540 
6540 
6550 
6560 
6570 
6580 
5590 
5500 
6610 
5620 
6630 
6640 
6650 
5659 
5570 
5680 
5570 
6700 
6710 
6720 
6730 
6741 
6750 
6760 
6770 


| „60 5ЕТ ۶لا‎ ٢.۶ USING D s PRESEA Е 
PRINT CH2$1(12) 
IF Tratio=A THEN CALL Plot2(E,D(*) ,x) 
CALL Distribution (D(#) „Јгатдо,Е) 
NEXT Tratio 
CN KEY 2 LABEL "HARD COPY" GOTO Hard 


ON KEY 3 LABEL "HENUE ' GOTO Seberit 
PRINT USING */////" 
PRINT ' YOU САМ GET A HARD COPY BY PRESSING KEY # 2° 
PRINT ' OR EXIT THIS ROUTINE KY PRESSING KEY $ 3" 
Standby:G0TO Standby 
Hard:X=1 IA VARIABLE TO CONTROL PLOTTER 
GOTO Begin 
Subexit: GINIT 


GRAPHICS OFF 
SUBEND 
SUB Shift 
COM /Readata/ B,P,H,Q5$[201,84$1201,7:$[201, Y1(x) BUFFER,Y2(x) BUFFER 
БИРО ЫЕ) 
PRINT CH8$(12) 
PRINT USING "///////" 
PRINT " RAW DATA IS BEING SHIFTED TO CORRECT FOR SMALL GAPS RETRKEEN SCANS" 
PRINT '" 
PRINT ° BE PATIENT, IT'S A LONG SET OF DO-LOOPS" 
| THERE ARE SOME SMALL GAPS BETWEEN SCANS AND THIS SUBROUTINE 
! SHIFTS THE DATA SQ THAT THE FIRST DIODE DATA POINT IS MOVED 
l TO THE VERY BEGINNING OF ITS 1424 BLOCK IN THE OVERALL ARRAY 
! .THE FIRST SET IS RIGHT ON, THE NEXT IS ONE OFF,THE THIRD IS 
| TWO OFF ,,50 FORTH. THIS MAY NOT MATTER WITH OUR RESOLUTION 
! AND IS PRORAPLY DUE TO THE MEMORY CARD CYCLING AT THE END OF 
| SOME SCANS. SEE THE CIRCUIT TIMING DIAGRAM IN THE THESIS. 
SELE 
CASE 4 IEXHAUST DATA HAS 4 SCANS 
М=3 
CASE 8 IMOTOR DATA HAS 8 SCANS 
M-7 
END SELECT 
FOR K=1 10 2 
IF K=1 THEN MAT E= 71 IOPPERATES ON NO-PARTICLE AND PARTICLE SETS 
IF К=а THEN MAT E= Y2 
FOR J-0 TO M 
IF 4-7 AND J{=3 THEN 6800  !One 4096 Block Doesn't Meed Shifting 
FOR 1=(J)*1024+1 TO (J+1) #1024 !ELOCKS OF 1024 
IF M=3 THEN L=ItJ lArray B Has the worst 
Problem With Shifting Data 
IF M=7 THEN L=I+1 Array D 1s Allways off by one 
IF L)(J+1)x1024 THEN L=(J+1) #1024 
!JUST TO AVOID PROGRAMMING ERROR AT THE END OF THE ARRAY 
E(T)=E(L) ! THIS SHIFTS THE DATA 
۱ IDEPENDING ON 'L',ARRIVED АТ 
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$780 
5790 
6800 
6810 
5820 
5830 
5840 
5850 
6860 
5870 
5880 
6890 
6900 
6910 
5920 
6730 
5940 
5950 
5960 
6270 
5989 
6790 
7009 
7010 
7020 
7030 
7049 
7050 
7060 
7070 
7090 
7090 
7100 
7110 
7120 


7130. 


7140 
7150 
7160 
7170 
7180 
7190 
7200 
7010 
72 

7230 
7240 
7250 
7260 
7270 


IBY OBSERVING RAW DATA 
NEXT I 
HEXT J 
IF К=1 THEN МАТ 
IF Kz2 THEN MAT 
NEXT k 
SUREND 
SUB Store 


ПЕ Е 
Ye= E 


THE LARGE ARRAY CONTAINING MULTIPLE SCANS HAS PEEN REDUCED 


COH /Hrdgaus/ Av2(*) 
DIM Data(1:1024) RUFFER 


TO A MEAN SCATTERING PROFILE BY AVERAGING AND FILTERING. IF 
YOU FEEL THAT STORING THIS DATA IS NECESSARY THIS ROUTINE 
DOES IT, IT SAVES DISK SPACE TO STORE THE RESULTS THER 
ELIMINATE THE RAW DATA IF YOU FEEL CONFIDENT THAT THE 
| REDUCTION IS THE BEST THAT CAN KE. IN OTHER WORDS, 
DO NOT PURGE A RAW DATA FILE UNLESS YOU ARE ABSOLUTELY SURE 
YOU WON’T WANT TO REDUCE IT AGAIN. 

lAve(*) IS THE REDUCED DATA 
!BUFFER USED TO TRANSFER TO DISK 


= = = o— оле тыз ---- °. зз 


MAT Data- Ave 
PRINT CHR$(12) 
PRINT USING *'///7//" 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


INPUT 44 
CREATE BDAT A$,512,16 

ASSIGN BDisk TO A$ 

ASSIGN BBuff TO BUFFER Data(¥) 


CONTROL 8Buf f ,3,1,8192, 


А SUGGESTED METHOD OF NAMING REDUCED DATA FILES IS AS FOLLOWS" 
Qo - 7 


EXHAUST ۰۳۰‏ سو 
Carros CALIBRATION , IF NO ’C’ THEN AN ACTUAL FIRING IS ASSUMED"‏ 
GLASS BEADS '‏ = 


IF NO ’C’ THEN ’G’ STANDS FOR ‘GAP’ PROPELLANT” 
A-—----RLUMINUM OXIDES 
PPP-----CHAMBER PRESSURE FOR RUN OR QTHER.,.. (HOZZLE 1۶۲۳ 
PARTICLE SIZE FOR CALIBRATION’ 

MMDD----MONTH, DAY OF RUN OR CALIBRATION” 

EXAMPLE: — MCA125JN12" 
MOTOR REAM CALIBRATION USING ALUMINUM OXIDE FROM 1 TO 25 MICRONS' 
ON JUNE 12" 

EXAMPLE: — XG550AU10 
EXHAUST BEAM DATA FOR GAP PROPELLANT AT 550 psi Oh AUGUST 10" 

PLACE A DISK IM THE RIGHT HAND DRIVE" 

ENTER THE FILENAME YOU WISH TO USE FOR THIS REDUCED DATA" 


! RUFFER IS FULL 


TRANSFER BRuff TO wDisk;COUNT 8192 
WAIT FOR EOT 8Disk 
ASSIGN 8Buff TU * 
ASSIGN #Disk TO # 


SUREND 
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7280 
7290 
7300 
731) 
7529 
73% 
7340 
7350 
7350 
7371 
7381 
7391 
7400 
7410 
7420 
7430 
7440 
7450 
7450 
7470 
7480 
7490 
7508 
7510 
7508 
7530 
7340 
7530 
7560 
7570 
7580 
7598 
7500 
7518 
7520 
7638 
7640 
7650 
7560 
7670 
7680 
7690 
7700 


SUB R 


eyview(Avl(*)) 


DIM Data(1:1024) BUFFER 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


* EACH REDUCED FILE CONTAINS ONE SCAN ONLY, YOU KUST KNOR IF IT" 

| IS EXHAUST , MOTOR CAVITY , ОК CALIBRATION DATA." 

* THE DISK WITH THE REDUCED FILE SHOULD BE IN THE RIGHT-HAND DRIVE” 
і ENTER THE FILENAME OF THE REDUCED DATA” 

қ TO BE REVIEWED” 

A$ 


ASSIGN Disk TO A$ 
ASSIGN 880۶۶ 10 RUFFER Data(*) 
CONTROL BDisk,5;1 


CONTROL 8Byff,3:1,0,1 
TRANSFER &Disk TO BBuff;COUNT 8192 


ITHIS IS AN EMPTY BUFFER 
!102418zNUMZER ОР FEYTES 


WAIT FOR EOT 8Disk 


ASSIG 


N @Disk TO * 


ASSIGN Buff TO + 


MAT À 


vi= Data 


SUBEND 


SUB F 
ШІ 2 
PRINT 
PRINT 
PRINT 


ilter (A(X), Fil) 
(1:1024) 

CHRS (12) 

USING *////////* 

' FILTERING THE SCATTERING PROFILE' 


FOR J=1 TO Fil 


FOR I=1 TO 1014 


B=] +] ITHIS IS A Symetric Moving Average 
С=1+0 ITYPE OF DIGITAL FILTER. EACH 
р=1+3 IDATA POINT IS EQUALLY WEIGHTED 
Р-1+4 IIN THIS CASE RUT THIS CAN RE 
6=145 ICHANGED IF ONE DETERMINES THAT 
H=I+6 IFEWER POINTS WITH UNEQUAL WEIGHTS 
K=I+? IWOULD RE FASTER OR GIVE BETTER 
[=]+8 IRESULTS, THIS TYPE OF FILTER WAS 
N=I+9 IUSED SINCE IT INTRODUCES NO PHASE 
P=] +10 ILAG (Anauiar Resolution), 
E(G)=(A(I)+A(B)+A(C)+A(D)+A(F)+A(G)+A(H)+A(K)+A(L)+A(N)+A(P))/11 
NEXT I 
MAT A= E 
NEXT J 
SUBEND 
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